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The influence of delamination on the progressing damage path and initial failure load 
in composite laminates is investigated. Results are presented from a numerical and an 
experimental study of center-notched tensile-loaded coupons. The numerical study includes two 
approaches. The first approach considers only intralaminar (fiber breakage and matrix 
cracking) damage modes in calculating the progression of the damage path. In the second 
approach, the model is extended to consider the effect of interlaminar (delamination) damage 
modes in addition to the intralaminar damage modes. The intralaminar damage is modeled 
using progressive damage analysis (PDA) methodology implemented with the VUMAT 
subroutine in the ABAQUS finite element code. The interlaminar damage mode has been 
simulated using cohesive elements in ABAQUS. In the experimental study, 2-3 specimens each 
of two different stacking sequences of center-notched laminates are tensile loaded. The 
numerical results from the two different modeling approaches are compared with each other 
and the experimentally observed results for both laminate types. The comparisons reveal that 
the second modeling approach, where the delamination damage mode is included together with 
the intralaminar damage modes, better simulates the experimentally observed damage modes 
and damage paths, which were characterized by splitting failures perpendicular to the notch 
tips in one or more layers. Additionally, the inclusion of the delamination mode resulted in a 
better prediction of the loads at which the failure took place, which were higher than those 
predicted by the first modeling approach which did not include delaminations. 

I. Introduction 

As the motivation increases to use composite materials in aerospace structures, researchers are 
investigating the sensitivity of these materials to damage. There are a number of sources of airframe initial damage. 
Two that have received considerable attention are delaminations and notches. A delamination can occur during 
manufacturing or from an impact on the structure. Also, delaminations can occur due to interlaminar stresses [1]. 
These stresses are generated due to mismatches in Poisson’s ratio and in-plane stress gradients. A through notch can 
result from a forklift or a rotorburst puncture of a fuselage, for example. Designers want to know whether damage 
will propagate under loading in order to assess continued airworthiness. Propagating damage can be intralaminar in 
nature (fiber breakage or matrix cracking within a layer of the laminate) or interlaminar in nature (delamination or 
disbonding of adjacent layers of the laminate). In many cases, the accumulation of damage involves both the intra 
and interlaminar modes. In such cases it is important to consider both damage modes in a progressive damage 
analysis (PDA) of the structure in order to accurately predict the damage path and failure load. 
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There have been numerous studies considering the effects of interlaminar damage in composites. Some 
researchers in the past have shown that stacking sequence has an influence on the initiation of delaminations [2-5]. 
Other researchers have studied and concluded that delaminations can have a significant effect on a laminate’s 
compressive strength. Progressive delamination analysis has been performed in Reference [6], to predict 
delamination growth between a composite skin and a stiffener using the Virtual Crack Closure Technique (VCCT) 
implemented within ABAQUS. In References [7-8] interface elements were used to model inter-laminar damage 
initiation and its progression, whereas in Reference [9], inter-laminar damage was modeled using a fracture energy 
based criterion for initiation and propagation of damage between two lamina. The authors in Reference [10] have 
utilized a cohesive element in the ABAQUS finite element code to predict the delamination response in a notched 
tensile-loaded laminated coupon. Intralaminar damage was not considered and no comparison with experimentally 
observed failure loads was mentioned. There have also been numerous studies considering the effects of 
intralaminar damage in composites. Intralaminar failure modes have been modeled in References [9, 11-14] using 
Hashin’s failure criteria [15]. Hashin’s failure criteria have been modified to account for tensile fiber failure due to 
transverse shear loading in References [16]. Qing, Chang and Starnes [17] used an intralaminar progressive failure 
algorithm in conjunction with ABAQUS to predict the damage patterns and tensile failure loads of notched 
laminates with co-cured reinforcing strips. Analytical failure loads agreed well with test data for a variety of notch 
sizes and reinforcing strip thicknesses when the damage mode was dominated by self-similar crack progression. 
When thicker reinforcing strips were used, splitting accompanied by delamination of the strip occurred when the 
damage reached the strip. This led to a 12% over-prediction in failure load. In the authors’ previous research [18] it 
was observed that a modeling approach that included only the intralaminar damage modes, underpredicted the 
failure load by 20% in an un-stiffened notched tension loaded laminate, where significant delamination was 
observed experimentally. These results support the inclusion of interlaminar effects in future code development. 
There have not been many studies where both intra and interlaminar damage modes are considered together. The 
modeling in Reference [7] considers the interaction between intralaminar matrix cracking and delamination by 
providing an empirical coupling relationship between matrix cracking and delamination. 

A majority of the progressive damage analyses have considered either interlaminar failure alone or 
intralaminar failures only and have not accounted for the interaction between the two, resulting in considerable 
discrepancies between the predicted and experimental results. Motivated by the need to more accurately predict the 
damage modes, and the damage progression in laminated structures, in the current work a PDA approach has been 
implemented that extends our previous work [18] by including both intralaminar and interlaminar damage modes in 
the prediction of damage progression. The methodology should be applicable to general damaged or undamaged 
laminated structures. The focus in this paper will be the application of the methodology to tensile loaded notched 
laminated coupons. 

In this paper, results are presented from a numerical and experimental study of center-notched tensile 
loaded specimens. The numerical study includes two approaches. In the first approach, the author’s previous PDA 
approach is utilized. It only considers intralaminar damage modes in calculating the progression of the damage path. 
Damage initiation and post initiation degradation of the properties of the failed elements associated with 
intralaminar failure modes have been implemented through the VUMAT subroutine in the ABAQUS finite element 
code. In the second approach, the effect of the interlaminar damage mode in addition to the intralaminar modes are 
considered in calculating the progression of the damage path. The interlaminar damage mode (delamination) 
predictions have been incorporated by using COHESIVE elements within ABAQUS. The numerical models 
represented the specimens used in the experimental study which consisted of 2-3 specimens each of two different 
laminates. Both laminates types consist of 10 plys with the same percentages of 0°, 45° and 90° plies but had 
different stacking sequences. The numerical results from the two different modeling approaches, which include 
load-end displacement, notch tip strains, and individual ply damage modes, are compared with each other and the 
experimentally observed results for both laminate types. In the second modeling approach, where delaminations are 
considered, the calculated delaminations are compared with those observed experimentally for both laminate types. 
An assessment of the comparisons reveals that the second modeling approach, where delamination damage modes 
are included, better matches the experimentally observed damage paths, which were predominantly splitting 
perpendicular to the notch tips in a number of layers, and the loads at which the splitting took place, which were 
higher than those predicted by the first modeling approach which did not include delaminations. 

II. Progressive Damage Analysis 

An analysis methodology for predicting damage initiation and propagation, and failure of laminated 
composite panels with a center notch is developed using a progressive damage analysis methodology implemented 
within the ABAQUS/explicit finite element analysis code. The use of an explicit analysis code for simulating the 
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progressive damage of statically loaded specimens is documented in the author’s previous work [18] but has much 
to do with the desire to predict fast running cracks in the future. The progressive damage analysis methodology 
includes intralaminar and interlaminar damage initiation and progression. Both geometric and material 
nonlinearities are included in the analysis. Material nonlinear behavior is associated with degradation of the 
material properties of the plies and degradation of the material in-between the plies to simulate the intralaminar and 
interlaminar damage mechanisms, respectively. Details of the methods used to simulate the intralaminar and 
interlaminar damage progression are presented in the next two subsections. 


A. Intralaminar Damage 

The intralaminar damage modes considered in the current study are: 

a. Fiber failure in tension and compression, and 

b. matrix failure in tension, compression and shear. 

The Hashin-Rotem unidirectional criteria [16] are used in the present study to predict intralaminar damage initiation. 
The ply stresses in the principal material directions are used in the unidirectional failure criteria. The failure criteria 
are expressed in terms of the in-plane stresses o ip and the strengths X, and Y and the shear allowable S. The in-plane 
stresses are cr n , cr 22 , and <Jn, where direction 1 is the fiber direction and direction 2 is the transverse direction in the 
ply. The strengths parallel and transverse to the fiber direction are denoted as X and Y, respectively, and the in-plane 
shear strength is *Sj 2 . The subscripts 'c' and ’f denote compression and tension, respectively. In the current study, in- 
situ strengths were used for Y t , Y c and S 12 [19]. 

The fiber- failure index, e f , in tension and in compression is defined as 
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The progressive damage analysis in ABAQUS/Explicit is implemented using a user- written subroutine 
VUMAT [20]. The VUMAT subroutine is invoked at each material point of an element to evaluate the failure 
criteria. When damage initiation is predicted at the material point, stresses are set to zero according to the 
degradation model presented in Table 1 to represent the degradation of the material associated with different damage 
modes. When the fiber failure criteria is satisfied, or, when any of the transverse or shear strains reach 0.5 at the 
material point, then that material point in the element is flagged for deletion. Once a material point has been 
deleted, it cannot be reactivated. The ABAQUS/Explicit code checks for elements where all of the material points 
have been flagged as deleted and removes these elements from further computations. 

B. Interlaminar Damage (Delamination) 

Interlaminar damage is simulated using the ABAQUS cohesive finite element COH3D8 at the ply interfaces. 
The interlaminar damage model is based on a continuum-damage-mechanics formulation of a stress-strain 
constitutive law for a thin resin layer that exists between adjacent composite lamina. The on-set of delamination is 
determined based on an interlaminar maximum stress obtained from a multi-axial stress criterion and the 
delamination growth is based on a critical fracture energy criterion. Damage is modeled as irreversible by including 
a damage parameter. The damage parameter is directly related to the dissipated fracture energy. The guidelines 
presented in Reference [21] have been used in the current study to calculate cohesive layer parameters such as mesh 
size, stiffnesses and the interfacial strengths of the interface (cohesive layer) between the plies. The equation 
presented in Reference [21] to define the stiffness of the cohesive layer in the mode-I direction is: 
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Where E 3 is the Young’s modulus of the laminate in the thickness direction, t is the larger of the sublaminate 
thicknesses above or below the cohesive layer and (X is a parameter that is much larger than 1 . The authors in 
Reference [21] recommend the use of a value equal to or greater than 50 for CC to obtain a stiffness of the cohesive 
layer which is small enough to avoid numerical problems such as spurious oscillations of the tractions in an element 
and also large enough to prevent the laminate from being too compliant in the thickness direction. Hence in this 
study a value of 50 has been considered. In calculating the stiffness, K n and K IU in the shear directions (mode-II and 
mode-III directions), E 3 is replaced by the shear moduli Gi 2 and Gi 3 of the laminate respectively. 

Interfacial strengths, required to simulate delamination, are determined based on the in-plane mesh size of the 
cohesive layer, the fracture toughness of the material and the number of elements in the cohesive zone using the 
following equation for mode-I failure: 
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where Uj is the interlaminar strength in the mode-I direction, G/ is the corresponding fracture toughness, N e is the 
desired number of elements in the cohesive zone and l e is the length of the finite element in the cohesive layer [21]. 
The other interfacial strengths, U n and U IU are computed in the mode-II and mode-III directions by replacing G/ in 
the above equation with the G n and G ni respectively. 

III. Experimental Description 

Two laminate stacking sequences, shown in Table 3 were considered in the progressive damage analysis 
study. Both laminates had the same percentages of 0°, 45° and 90° degree plies, but the stacking sequence was 
varied to study the evolution of delamination at different interfaces. The geometric description of the center-notched 
laminate specimens is shown in Figure 1. The specimens have a test section of 6 inches x 3 inches. The specimens 
have a 2.75 inch long tab section on each end, which fits into the grips of a servo-hydraulic testing machine and the 
specimens have a notch dimension of 0.75 inches x 0.09375 inches. 

Four axial strain gauges were mounted on the specimens to measure strain data during the loading. The 
strain gage locations on the specimens are shown in Figure 2. EA-06-070LC-350 type strain gages were mounted 
on the notch tips. CEA-06-187UW-350 type gages were mounted at the other locations. Two extensometers of 4 
inches span length were installed on the two vertical longitudinal edges of the specimens, with the ends equidistant 
from the horizontal notch line to provide a direct measurement of the elongation throughout the loading. Figure 3 
shows one of the specimens mounted in the test machine. The specimens were statically loaded in tension at room 
temperature, under displacement control at a rate of 0.005 inches/minute, and the test data were acquired every half 
second. 

The specimens are fabricated from T800/3 900-2 pre-preg tape, a toughened epoxy CFRP composite 
material. Stiffness properties and strength allowables are given in Table 2. The laminate sequences and their 
corresponding thicknesses are defined in Table 3. Three specimens of laminate- 1 and two specimens of laminate-2 
were loaded above the initial failure load but lower than the ultimate failure load of the specimens. The observed 
initial failure load was splitting of some of the plies at the notch tips, with splits propagating perpendicular to the 
notch direction all the way to the grips. The intralaminar splitting was accompanied by delaminations. The ultimate 
failure load, characterized by the separation of the specimen into two or more pieces, occurred at a load a few 
thousand pounds higher. Only one specimen of laminate-2 was loaded to ultimate failure, at which point, post test 
investigation of the specimen was meaningless. 

IV. Finite Element Model Description 

Two finite element models, one a 2-D model and the other a 3-D model were developed to perform the 
PDA using the ABAQUS/Explicit code. The 2-D, 4-node shell model, used in the predictions of intralaminar 
damage modes only, is shown in Figure 4. The finite element model includes only the test section, which is of 3.0 
inches wide and 6.0 inches long with a center notch of 0.75 inches x 0.09375 inches. The 2-D finite element model 
of the specimen is discretized using 4-node quadrilateral elements whose size is approximately 0.08 inches x 0.08 
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inches. The finite element type used in the intralaminar progressive failure analyses is the ABAQUS S4R and the 
model consists of 3012 nodes and 2888 elements for each laminate sequence. 

A 3-D model was developed to simulate coupled inter- and intralaminar damage modes. The thickness of 
this model is divided in to 5 sections with a cohesive layer between each of the sections. The pictorial 
representation of the thickness discretization of both of the laminates is shown in Figure 5. In case of laminate- 1 , 
the 1 st and 5 th sections consist of [45°] plies, the 2 nd and 4 th sections consist of [07-4570°] plies and the 3 rd section 
consists of [90790°] plies. The cohesive layers were included at the 0790° and the 0745° interfaces where 
delamination is most likely to occur due to Poision’s ratio mismatch and the coefficient of mutual influence of the 
adjacent plies respectively. In the case of laminate-2, the 1 st and 5 th sections consist of [070°] plies, the 2 nd and 4 th 
sections consist of [45°] plies and the 3 rd section consists of [-457907907-45°] plies. The cohesive layers were 
modeled between 0° and 45° plies and also between 45° and -45° plies, interfaces which were unique to laminate-2. 
The stacking sequences for the laminates including the cohesive layers are presented in Table 4. Each section of the 
3-D model is discretized using one ABAQUS SC8R (8-node shell continuum) element in the thickness direction. 
The assemblage of the five sections and four cohesive layers through the thickness defines the laminates. SC8R 
elements of size 0.08 inches x 0.08 inches, were adopted to discretize the in-plane direction. The cohesive layers 
were modeled using a zero thickness ABAQUS 8 -node cohesive element known as COH3D8. An in-plane 
discretization of 0.02 inches x 0.02 inches, was maintained in the cohesive layers which require a finer mesh. The 
interfacial boundaries between the sub-laminate and the cohesive layers were connected using “Tie” boundary 
conditions to accommodate the different mesh sizes. 

In the 2-D model, one of the transverse (horizontal) edges is fixed in all the translational and rotational 
degrees of freedom about the x, y and z directions. The other transverse edge is fixed in the y and z translational and 
x, y and z rotational degrees of freedom but a uniform end-displacement was applied in the x direction. In case of 
the 3-D model, one of the transverse (horizontal) edges is fixed in all the translational degrees of freedom (x, y and z 
directions). The other transverse edge is fixed in the y and z translational degree of freedom but a uniform end- 
displacement was applied in the x direction. In both the 2-D and 3-D models, the two long (vertical) edges are left 
free. 

The 3-D model requires interfacial properties of the cohesive layer to simulate delaminations. Hence, the 
properties such as interlaminar stiffnesses and strengths are evaluated using equations 3 and 4. In equation 3 the 
thickness t was determined based on the maximum number of plies between any cohesive layers. Hence the values 
of t equal to 3 ply thicknesses and 4 ply thicknesses were considered for laminate- 1 and laminate-2 respectively, to 
calculate interfacial stiffnesses. The interfacial strengths were calculated by adopting two finite elements in the 
cohesive zone [21] with the length of the element being 0.02 inches. The interlaminar strengths and stiffnesses 
values are presented in Table 5 for both laminate- 1 and laminate-2. 

The progressive damage analysis was performed by applying a uniform end displacement of 0.05 inches 
and 0.06 inches along the 0° (fiber) direction for the 2-D and 3-D models respectively, over a period of 0.01 
seconds. The end displacement was applied in a ramp fashion and a mass density of 0.000154 lbs-sec 2 /in 4 for the 
composite material was used in the analysis. During the entire simulation, the kinetic energy of the model was 
checked to ascertain that it did not exceed 10% of the total strain energy of the model in accordance with the 
ABAQUS user’s manual guidelines for simulating static problems using the explicit finite element code [20]. This 
was necessary to ensure that while simulating a pseudo-static analysis, significant dynamic effects were not 
introduced while using the fictitiously small total loading time required for computational efficiency in an explicit 
analysis. 


V. Results and Discussion 

Results including the load-end displacement curve and the notch-tip strain, obtained from both the 
experiment and the 2-D and 3-D progressive damage analyses are presented. C-scan pictures of the damaged 
specimens and the simulated interlaminar and intralaminar damage paths are also presented in Figures 6-17. The 
axial stress distribution, before and after delamination initiation in the 0° ply, obtained from the 3-D PDA is 
presented for laminate-2. 

A. Laminate-1 

Three specimens of this laminate sequence were loaded to 15,300.0 lbs in tension. In Figure 6, the load- 
end displacement curve obtained from the experiment, and the curves predicted from the 2-D and 3-D PDA have 
been presented. In this figure, it can be seen that the 2-D PDA predicts complete failure of the specimen at 13,613.0 
lbs which is represented by the red curve in the figure, indicating the loss of load carrying capability of the 
specimen. In contrast, the curve predicted from the 3-D PDA, shown in green, indicates a small change in the slope 
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at 13,330.0 lbs., but shows that the laminate continued to carry load, which is more consistent with the experiment, 
shown in blue. Also, at 13,300.01bs, it was observed that the 3-D PDA predicts matrix damage and delaminations 
around the notch tip, and fiber damage in the -45° plies. At 15,316.0 lbs the simulation was stopped as was the case 
with the experiment. The load vs. end displacement prediction from 3-D PDA and the experimental data agree 
reasonably well up to the maximum applied load. In Figure 7, the strain predictions from 2-D and 3-D PDA, and the 
measured strain at the notch-tip are presented. It can be observed in this figure that the measured strain remains 
linear only within a short load range and then becomes nonlinear as the specimen is loaded further. The high strain 
value, which is in excess of 3%, is well above the fiber failure strain. This implies that this is a surface phenomenon 
at the notch-tips and needs additional study. 

In Figure 8, the damage mode predictions from the 2-D PDA are presented. In this figure all of the plies 
except the 90° ply in the laminate failed in the fiber damage mode. However the 90 degree plies failed in the matrix 
damage mode. Hence, it is obvious from these pictures that the 2-D PDA predicts a through crack, initiating at the 
notch tips with propagation in a self-similar fashion all the way to the free edge when all of the damage modes are 
integrated through the thickness. The red and blue colors in the picture indicate the damaged and the undamaged 
areas in the plies respectively. In Figure 9 the intralaminar damage predictions from the 3-D PDA are presented. 
The dominant damage mode in all of the plies except the -45° plies is in the matrix damage mode with splitting 
occurring perpendicular to the notch tip. The -45° plies in the laminate failed in the fiber damage mode which in turn 
caused a noticeable change in the stiffness of the laminate, which also can be observed in Figure 6. However the 
load carrying plies (0° plies) did not show fiber damage when the simulation was stopped indicating additional 
unused load carrying capacity. In Figure 10 the interlaminar damage (delamination) predictions at the 4570° 
interface and the 0790° interface are presented. The red and blue/green colors in the picture indicate complete 
delamination and no delamination respectively. Comparing the simulated delamination with the one recorded by the 
C-scan in Figure 11, it is clear that there is a good qualitative agreement between them. Also in Figure lib, the 
popping of the 45 degree ply starting from the notch-tip, splitting in a “Z” fashion due to delamination at the 45/0 
interfaces is presented. 

B. Laminate-2 

Three specimens of this laminate sequence were loaded in tension. Two specimens were loaded to 
21,280.0 lbs and the other one was loaded to 23,820.0 lbs which caused the specimen to fail completely. The 
experimental results presented in Figures 12-17 correspond to the specimens which were loaded to 21,280.0 lbs. In 
Figure 12 the load-end displacement curves obtained from the experiment, 2-D and 3-D PDAs are presented. The 
curve predicted by the 2-D PDA indicates a complete self-similar failure at 13,618.0 lbs. However, the curve 
predicted by the 3-D PDA indicates a small change in the slope from the measured load-end displacement curve and 
further load carrying capacity to 14,387.0 lbs. Unfortunately, at this load level the analysis terminated due to an 
excessive distortion in one of the elements. Therefore it was not possible to see the additional load carrying capacity 
observed in the experiment. In Figure 13, the axial strain predictions at the notch-tip from the 2-D and 3-D PDA’s 
and the experimental results are presented. The important observation is that the strain values predicted by the 3-D 
PDA are matching very well with the experimental data, which was not the case for laminate- 1 . And moreover, as 
laminate-2 has 0° outer plies, the nonlinearity in the measured load-strain curve at the notch-tip is significantly 
reduced as compared to the one shown for laminate- 1, which has 45° outer plies. From this observation, the authors 
feel that the non-linearity of the strain at the notch-tips observed in laminate- 1 is a local phenomenon. In Figure 14, 
the axial stress distribution in one of the 0° plies (load carrying plies) before and after the delamination at the 0745° 
interface is presented. It can be clearly noticed from this figure that the delamination decreases the notch-tips stress 
and hence increases the load carrying capability of the specimen by redistributing load in the load carrying 0° plies. 

In Figure 15, the damage mode prediction from 2-D PDA is presented. In this figure all of the plies except 
the 90° ply in the laminate failed in the fiber damage mode again in a self-similar manner. However the 90° plies 
failed in the matrix damage mode. In Figure 16, only the matrix damage predictions from the 3-D PDA are 
presented as there was no fiber damage present in any of the plies. In Figures 17 and 18 the delaminations predicted 
at the 0745° and 457-45° interfaces from the 3-D PDA and the delamination observed experimentally from the C- 
Scan are presented. It can be observed that the predicted delamination shape is in reasonable agreement with the 
measured one. However, hardly any delamination was predicted at the 457-45° interface. 

VI. Conclusions 

The influence of delamination in altering the damage modes, the damage path and the initial failure load of 
notched composite specimens has been successfully demonstrated. It is evident that a 2-D PDA with intralaminar 
failures only, predicts complete failure of the specimen at a lower load and with self similar cracks emanating from 
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the notch tips for laminate- 1 and laminate-2. The interlaminar damage mode (delamination) has been predicted for 
two tensile loaded notched laminate sequences using the 3-D PDA. From the 3-D PDA, delaminations were 
predicted between plies, which in turn increased the load carrying capacity of the specimen which is evident in 
Figures 6 and 12. The authors’ believe that the predicted and experimentally observed increase in load carrying 
capacity is due to load redistribution around the notch tips, shown in Figure 14, caused by energy dissipation during 
the delamination process. The intralaminar damage mode and damage path predictions from the 2-D and 3-D PDA 
were also different in most of the plies for both the laminates. Hence, these results suggest that the delamination 
alters the damage modes at the notch-tips considerably and without modeling delamination, one would arrive at the 
wrong damage modes and ultimately the wrong initial failure load for laminates prone to delamination. 
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Table 1. Degradation model for ABAQUS/Explicit VUMAT. 


Damage mode 

Stresses set to zero 

Fiber failure in tension and compression 
Matrix cracking in tension, compression and shear 

All stresses 
^22, <*12 



Table 2. Material properties and strength data. 

Properties 

T800/3900-2 
(X 10 6 psi) 

Description 

Eii 

23.2 

Young’s modulus in fiber direction 

E22 

1.3 

Young modulus in the transverse direction 

G 12 

0.9 

In-plane shear modulus 

x, 

0.412 

Tensile strength in the fiber direction 

x c 

0.225 

Compressive strength in the fiber direction 

Y, 

0.00872 

Tensile strength in the transverse direction 

Y c 

0.0243 

Compressive strength in the transverse direction 

S 12 

0.01376 

In-plane shear strength 

V12 

0.28 

Poission’s Ratio 

Gi 

0.86 in-lbs/in 2 

Energy release rate in mode I 

Gu, Gm 

2.7 in- lbs/in 2 

Energy release rate in mode II and mode III 


Table 3. Specimen Specifications of 4-Node shell model. 

ID 

Ply Percentages 
[0’s/45’s/90’s] 

Laminate Stacking Sequence 

Thickness 

(inches) 

Laminate - 1 
Laminate - 2 

[40%/40%/20%] 

[40%/40%/20%] 

[45/0/-45/0/90] s 

[0 2 /45/-45/90] s 

0.065 

0.065 

Table 4. Specimen Specifications of 8-Node shell model. 

ID 

Ply Percentages 
[0’s/45’s/90’s] 

Laminate Stacking Sequence 

Thickness 

(inches) 

Laminate - 1 
Laminate - 2 

[40%/40%/20%] 

[40%/40%/20%] 

[45/decoh/0/-45/0/decoh/90] s 

[0 2 /decoh/45/decoh/-45/90] s 

0.065 

0.065 


Table 5. Cohesive Layer Parameters. 


Parameters 

Cohesive layer 

Description 

U, 

4969 psi 

Interlaminar shear strength in model I fracture mode 

u„, U m 

8805 psi 

Interlaminar shear strength in mode II and mode III fracture mode 

Kj 

3.4E+9 

Laminate - 1 

Kn, K UI 

1.2E+9 

Laminate - 1 

Kj 

2.5E+9 

Laminate - 2 

Kn, K m 

8.8E+8 

Laminate - 2 
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Figure 1. Center-notched Tensile Specimen. 
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Figure 2. Strain gauge locations on the specimen. 
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Figure 3. Test Setup. 


U 



R x =R y =R x =0 

Figure 4. Finite-Element Mesh. 
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Cohesive 

elements 




Cohesive 

elements 


Laminate - 2 


Figure 5. Laminates showing the location of cohesive layers. 



Figure 6. Load vs. End-Displacement of Laminate-1. 
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Figure 7. Load vs. strain in Laminate-1. 



0 Degree ply 45 Degree ply -45 Degree ply 90 Degree ply 

Fiber Failure Fiber Failure Fiber Failure Matrix Failure 

Figure 8. Intralaminar damages obtained from 2-D PDA of Laminate-1 at P x = 13613.0 lbs, 



0 Degree ply 45 Degree ply -45 Degree ply 90 Degree ply 

Matrix Failure Matrix Failure Fiber Failure Matrix Failure 

Figure 9. Intralaminar damages obtained from 3-D PDA of Laminate-1 at P x = 15316.0 lbs. 



a) 45/0 Interface b) 0/90 Interface 

Figure 10. Delamination predictions for laminate-1 at P x = 15316.0 lbs. 
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Figure 11. C-Scan picture showing delamination and a digital photo showing debond in laminate-1 at 

P x = 15300.0 lbs. 



Displacement, A x , of a 4-in. long section, inches 


Figure 12. Load vs. End-displacement of Laminate-2. 



Figure 13. Load vs. strain in Laminate-2. 
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190,100.0 psi 


167,200.0 psi 




a) P x = 8,371.0 lbs 


b) P x = 13,411.0 lbs 


Figure 14. Axial Stress distribution in 0° ply of Laminate-2 before and after delamination at 0/45 interface. 



0 Degree ply 45 Degree ply -45 Degree ply 90 Degree ply 

Fiber Failure Fiber Failure Fiber Failure Matrix Failure 


Figure 14. Intralaminar damages obtained from 2-D PDA of Laminate-2 at P x = 13618.0 lbs. 



0 Degree ply 45 Degree ply -45 Degree ply 90 Degree ply 

Matrix Failure Matrix Failure Matrix Failure Matrix Failure 


Figure 15. Intralaminar damages obtained from 3-D PDA of Laminate-2 at P x = 14387.0 lbs. 
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a) 0/45 Interface 


b) 45/-45 Interface 


Figure 16. Delamination predictions for laminate-2 at Px = 14387.0 lbs. 



Figure 17. C-Scan picture of the damage in specimen of Laminate-2 at P x = 21280.0 lbs. 
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